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Abstract 

The presence of fields with negative mass-squared typically leads to some form of instability 
in standard field theories. The observation that, at least in the light-cone gauge, strings prop- 
agating in plane wave spacetimes can have worldsheet scalars with such tachyon-like masses 
suggests that the supergravity background may itself be unstable. To address this issue, we 
perform a perturbative analysis around the type IIB vacuum plane wave, the solution which 
most obviously generates worldsheet scalars with negative mass-squared. We argue that this 
background is perturbatively stable. 
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1 Introduction 



It has been known for some time that plane-fronted, parallel gravitational waves (pp-waves) are 
exact solutions of string theory Jl], [|. The same is true of pp-wave backgrounds with a non- 
trivial dilaton and null Neveu-Schwarz-Neveu-Schwarz (NS-NS) 3-form field strength [||, [J. These 
backgrounds solve the /3-function equations to all orders in q'. A particular class of pp-wave is 
the exact plane wave and in recent months the subject of plane waves in string theory has become 
an intense area of research, the reasons for this renaissance being as follows: after the discovery 
of maximally supersymmetric plane waves in eleven-dimensional [Q, || and ten-dimensional type 
IIB H supergravity, it was observed 0, ^] that these solutions can be thought of as Penrose limits 
of the AdSp x S D ~ P vacua of the respective theories. Despite the fact that the ten-dimensional case 
involves the Ramond-Ramond (R-R) 5-form, it was then shown || that superstring theory on 
this background can be solved in the light-cone gauge. Finally, a sector of N = 4 super- Yang- 
Mills (SYM) which is dual to string theory on this background was identified by Berenstein et 
al p|, for the first time allowing for truly stringy tests of, at least a special limit of, the AdS/CFT 



correspondence [12, 13, 14]. 

The Penrose limit [15] of any solution of the ten-dimensional supergravity theories [16] gives 
a plane wave, which preserves at least one-half of the supersymmetries, but sometimes more: for 
example, the Penrose limit of the AdSs x T 1,1 geometry, dual to a certain J\[ = 1 superconformal 
gauge theory [17], gives the same maximally supersymmetric plane wave [18, 19, |20|. Since this 
latter work, Penrose limits have been taken of a whole host of geometries which arise in string 
theory. Of most relevance to us are the plane waves one obtains by taking Penrose limits of the 



geometries dual to non-conformal gauge theories: the Pilch- Warner geometry [21, 22, 23 1, dual to 
an M = 1 supersymmetric renormalization group (RG) flow between N = 4 SYM and an M = 1 



superconformal infra-red (IR) fixed point 24, [2^, |2£|]; the near-horizon limits of the Dp-brane 
geometries for p / 3 p5[| ; and various other relevant backgrounds ||, 20, 3C, 31, |3q| . 

Issues concerning holography, even in the maximally supersymmetric example, remain mysteri- 
ous [33, [3^, 36, 37]. Although both this |Q and more general plane waves |^7f[ , have a boundary 



which is a null line, it is not clear as to what the holographic theory actually is. It is thus even 
harder to address more general issues regarding RG flows in the plane wave limit, though in this 
context they have been discussed in (29], |33|, 



I |30|, 38, 31]. Certain general statements can be made, 

however: the Penrose limit of a geometry dual to a non-conformal gauge theory results in a plane 
wave with a time-dependent profile. In light-cone gauge string theory, it is precisely this profile 
which gives a mass to the worldsheet scalars so, in this case, one generically finds simple harmonic 
oscillator modes with time-dependent masses [^5j| . 

Moreover, the Penrose limit of the near-horizon Dp-brane geometry for p ^ 3 gives rise to plane 
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waves with profiles of the "wrong sign"; that is, which generate worldsheet scalars with negative 
(mass) 2 p5[| P|. One might think that such terms signal instabilities although, since the "mass" of 
the worldsheet scalars is really only an artifact of the light-cone gauge |25|, [HJ , this naive thought 
may not be correct. Either way, the issues involved warrant further investigation. 

Already in |2S| j, it was noted that, for a certain range of parameters, one can generate what 
appear to be stringy instabilities in a plane wave background with non-trivial 3-form and 5-form 
field strengths. In this case, however, all the mass terms have the correct sign and the zero modes of 
the string are unaffected; one would thus not expect to see classical instabilities of the background. 
In fact, backgrounds with a non-trivial 3-form will prove to be too complicated to analyse in any 
detail here. Even in the maximally supersymmetric case, however, the dimensionless string theory 
Hamiltonian is |Ll|, [k| 

H = Vn 2 + M 2 N n , (1.1) 

n 

where M denotes the dimensionless mass of the worldsheet scalars and N n is the occupation number 
for each mode at level n. It would appear that switching the sign of M 2 will give rise to instabilities 
not only for stringy modes with but also for the zero modes themselves^]. And it is precisely 

these latter potential instabilities which one might expect to observe in a standard classical analysis. 

Our purpose here is to assess the perturbative stability or otherwise of plane waves. That is, we 
will consider whether there exist solutions of the linearized field equations which grow exponentially 
in time. As far as we are aware, this is the first study of classical stability for non-static, though 
possibly stationary, spacetimes. Although we are only able to explicitly solve the relevant equations 
for the plane waves with constant profiles (the so-called Cahen-Wallach spaces |^(|), this class of 
spacetimes certainly includes cases for which the (mass) 2 of the worldsheet scalars is negative. We 
will show that this class of plane waves does not exhibit classical instabilities. It would also seem 
that scalar field theory in this background is well-defined, in the sense that it too is stable [O]. 

In the following section, we briefly review pp- and plane waves in an arbitrary number of 
dimensions, commenting on some general properties relating to the issue of stability, including 
supersymmetry and particle production. We study the string theory dispersion relation in the ten- 
dimensional vacuum plane wave in section 3, to assess in what sense it shows signs of instabilities. 
In section 4, we compute the linearized vacuum field equations of type IIB supergravity, specialising 
to the light-cone gauge to restrict ourselves to purely physical degrees of freedom. The resulting 
equations are analysed in section 5. Our analysis of the metric perturbations is applicable to plane 
waves in any dimension, and that of the dilaton and form-field perturbations generalize fairly 
readily. We conclude in section 6. 

1 Further examples of the occurence of modes of negative (mass) 2 in Penrose limits of other supergravity solutions 
are studied in @ 0. 

2 Of course, upon switching the s ign of M 2 , one loses maximal supersymmetry and it is not obvious that the 
Hamiltonian in this case is given by (1.1) with a simple sign difference. This is one of the issues we will address. 



3 



2 The spacetime 

Consider the vacuum pp-wave solution in an arbitrary number of dimensions (in the following 
section we will specialize to waves of type IIB supergravity). The metric 

ds 2 = 2dudv + H(u,x)du 2 + dx i dx i (i = 1, . . . , D - 2), (2.1) 

has curvature 



Ruiuj — ~7;didjH, R uu — — -didiH, (2.2) 



-didjH, R uu — — -( 
and scalar Laplacian 

□ = 28 u d v - H(u, x)d 2 v + m, (2.3) 

where we should note that we will not distinguish between upper and lower transverse indices. In 
general, the only isometry of this spacetime is that generated by the null Killing vector d/dv. 

Taking e = for null geodesies, e = 1 for timelike geodesies, and with A denoting the affine 
parameter, this isometry gives rise to a conserved energy, E: 

u(\) = E\ + uq. (2.4) 

The constraint L = — e is used to solve for v : with a dot denoting differentiation with respect to A, 
we have 

v(X) = J (e + H(u, x)E 2 + x'x 1 ) dX + v , (2.5) 
where the transverse coordinates satisfy 

x* = \diHE 2 . (2.6) 

For the plane wave, with H(u,x) = —Aij(u)x l x J , one can solve this equation for arbitrary initial 
conditions x l (0) and x l (0). 

The general pp-wave is singular, albeit in a special sense since all curvature invariants vanish 
identically |, g. To see this, one constructs an orthonormal frame which is parallelly propagated 
along a time-like geodesic, with tangent vector 

T a ={E,v,± i ). (2.7) 

The unit normals 

<) = fa, ~Ai^i > ™ a = ( E > * + 4 i* ), (2-8) 



all obey (T • Vn) a = 0, as required for a parallelly propagated frame of reference. The non-zero 
components of the Riemann curvature with respect to this frame of reference, 

E 2 

HT)(i)(T)U) = ^bcd-L n ({) i n {j) = - 



R (Tm(T)U) = RabcdT a n b {i) T c n d {j) = -— 8^11, (2.9) 



R{m)(i)(m)(j) = Rabcdm a n b {i) m c n d {j) = — —didjH, (2.10) 



E 2 
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will generically diverge. If such a divergence occurs at some point in the spacetime, then an observer 
will encounter infinite tidal forces at that point, and the space is geodesically incomplete^]. The 
only non-singular pp-wave 0, |3| is the plane wave, for which 

H(u,x) = -Aijiu^xl, (2.11) 

where Aij(u) is a smooth function of u. (For constant Aij, this is a Lorentzian symmetric space 
[10, H|.) In this case, 

R (T)(i)(T)(j) = R(m)(i)(m)(J) = E 2 Aij(u), (2.12) 

and the tidal force is attractive for positive eigenvalues of Aij and repulsive for negative eigenvalues. 
It might appear that this notion of a singularity is observer-dependent, but this is not the case since 
the tidal force, divergent or otherwise, would be experienced by all time-like observers. Indeed, 
one can show that the metric function Aij{u) encodes diffeomorphism-invariant properties of the 
spacetime [41]. 



The nature of the tidal force is seen in the geodesies also. The vacuum field equation is simply 

txA(u) = 0, (2.13) 

and, for constant Aij, we can always diagonalize, taking Aij = X^Sij. Then there must be at least 
one negative eigenvalue, and we can write = i^fj)- The geodesic equations for this simple case 



are easily solved: u and v are as in (2.4) and ([2.5]), and 

x* = A 1 cosGu (i) £A) + B i sin(u (i) £A) (A (i) = ^. } ), (2.14) 

x l = A i exp(^ (4) £A) + B l expi-fi^EX) (A (i) = (2.15) 

The attractive tidal force in the first case gives oscillatory behaviour around the origin, and the 
repulsive tidal force in the latter case gives geodesies which are pushed off to infinity. In this sense, 
the vacuum spacetime exhibits an instability (and we will see the same behaviour in the zero modes 
of the string in this background below). This behaviour is reminiscent of de Sitter space, in which 
particles also accelerate forever and, just as in that case, the repulsive tidal force in the directions 
associated with negative eigenvalues does not necessarily give rise to instabilities in the usual senseQ. 
In the case of de Sitter space, stability is guaranteed by the "cosmic no hair" conjecture [48, 49, 50], 
and one might wonder if there is a similar argument for the case of the plane wave. 

3 It is amusing to note that one encounters these "pp-curvature singularities" in the solutions describing pp- 
waves on fundamental strings [ fi^ , ff3| , [ii] , ^5[ , and also in non-dilatonic "brane-wave" solutions of ten- and eleven- 
dimensional supergravity |46f . In that case, even the exact plane waves give rise to such curvature singularities at the 
would-be horizon of the brane. Presumably, the related solutions considered in ^] have a similar global structure. 

4 The Penrose limit of pure de Sitter space is just flat space, as for pure anti-de Sitter space Although such 
a spacetime is not a solution of any known supergravity theory, if we consider the Penrose limit of de Sitter space 
times a sphere, along a geodesic which has angular momentum on the sphere, then we do find negative eigenvalues 
coming from the de Sitter directions. 
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Since all such plane waves are at least one-half supersymmetric, one might be tempted to 



argue that they must be stable. However, the extent to which the standard arguments pi], 52, 



53, |54| concerning stability and supersymmetry apply here is unclear. After all, the spacetime is 
not asymptotically flat. Moreover, the supersymmetries common to all such plane waves do not 
commute with the Hamiltonian. One can thus write down one-half supersymmetric configurations 
which nevertheless do not solve the relevant field equations However, demanding extra or, in 



the language of [56, 57], supernumerary, supersymmetries does restrict one to solutions of the field 
equations; and indeed it would seem that any such solution which preserves more than one-half of 
the supersymmetries does not give rise to the potential problems alluded to above. For example, 
it is fairly easy to see that waves with time-dependent profiles never have such supernumerary 
sup er symmetries . 

More explicitly, the integrability condition for the existence of Killing spinors is 

[A,, A]e = \didjH TjT+e = 0, (2.16) 

and multiplying this from the left with Tj gives 

dAH T+e = 0, (2.17) 

where didiH = is the field equation. One can thus have the standard supersymmetries, satisfying 
r + e = 0, without solving the field equation. Supernumerary supersymmetries, however, have 
r + e 7^ so, if any of these are to exist, the field equation must be satisfied. Of course, in this 
vacuum example, there are no supernumerary supersymmetries anyway. In the more general case 
with background fluxes, it is true that the existence of supernumerary supersymmetries implies the 
field equations, but this is harder to see. One has to consider the Killing spinor equation directly, 
in addition to the variation of any other fermionic fields, for example the dilatino in type IIB 
supergravity. An analogous study of the integrability condition for the existence of Killing spinors 
in this more general case is not enough^]. 

Further general comments are in order with respect to particle production in these backgrounds. 
It is easy to see that Gibbons' original argument |58|| — that there is no production of massless 
scalars in the sandwich wave — can be generalized to any pp-wave (see, e.g., [§, |3|). Given that 
d/dv generates an isometry, we can always decompose the modes of a scalar field in the pp-wave 
background with respect to v , as 

4>{u,v,x) = e ip - v 0(u,x). (2.18) 
To determine whether there is particle production, we need to consider whether the "in" modes 



3 We thank the anonymous referee for pointing out the problems with our original argument in this case. 
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in the asymptotic past are orthogonal to the "out" modes in the asymptotic future^. The Klein- 
Gordon scalar product is 

(0, = -* / d£ a ^ d a t/;*, (2-19) 
Jt, 

where S is the null surface u = 0, with unit normal d/dv. Despite the fact that these spacetimes 



are not globally hyperbolic (see also |p7[| ), this is a good substitute for a Cauchy surface [58]: 



it fails to capture only those geodesies moving parallel to the wave, with u = constant, and these 



will never interact with the wave itself. Explicitly, we have [58] 



(0,^) = -i / dvd D - 2 x (f)d v *p*, (2.20) 



D-2 

£ 

so the mode functions are always orthogonal, since they are simple exponentials in v. For any 
pp-wave, the null Killing vector d/dv allows for a preferred notion of frequency, so there can be no 
particle production and the resulting Bogoliubov coefficients must vanish. 

One might then argue that scalar fields in the pp-wave background are always stable, since 
presumably an instability in such a field will cause runaway particle production. What is less clear 
is whether this implies stability of the background itself (some comments in |]60| seem to be make 



this connection). One should also consider what effects interactions have [41]. 

3 String theory on the plane wave 

Of course, the zero modes of a string in the plane wave background will follow the same geodesies 
as above but we should also consider the non-zero modes to see what effect, if any, they have. 
Moreover, it would appear from (|l.l| ) that the zero modes themselves will give rise to instabilities 
for negative eigenvalues of the metric function Aij, and it is this issue that we want to address 
here. Ultimately, we are interested in whether or not there are tachyonic modes of the string (or 
particle) in the plane wave background. Since the notion of "mass" in this case is not the same as 
in flat space, we will ask the question as to whether there is faster-than-light propagation of any 
modes. It would appear that there is though, again, in a similar way as for de Sitter space. 

Since we are still working with the vacuum solution, the worldsheet fermions are unchanged 
from the usual flat space results. In the light-cone gauge, they are affected only by background 
fluxes ||10|| . Indeed, since the vacuum backgrounds are necessarily only one-half supersymmetric, 
there are no linearly realized worldsheet supersymmetries at all — after fixing the kappa-symmetry, 
only supernumerary supersymmetries are linearly realized |39]. Then it should not be surprising 
that, whereas the worldsheet scalars are massive, the fermions will still be massless. We will thus 
consider only bosonic string theory here. 



In the case of the sandwich wave, in the asymptotic past and future, the spacetime is flat and the obvious "in" 
and "out" states can be constructed |pq| . More generally, one can imagine |^5| an RG flow from some constant Aij to 
some other constant Aij in which case one can define "in" and "out" states in these stationary regions of spacetime. 
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The action 



S = ~4^J ^^^ haPd ^ a ^X b g ab , (3.1) 

with background metric 

ds 2 = 2dudv - ^-(u^Vdu 2 + di'dz* (i = 1, . . . ,8), (3.2) 
has the energy-momentum tensor 

T aP = ~ (d a X a dpX b g ab - l -h aP dX a ■ 8X b g ab ^ , (3.3) 



and canonical momenta 



P_ = |i = ij7, (3.4) 
dV it 



f)T 1 

P + = -^ = -(V-A ij X i Xiu), (3.5) 



Pi = ■S L = k* i - ( 3 - 6 ) 

dx* i 2 y ' 

where I 2 = 2na' . The V equation implies that DU = 0, so taking U = U(r), we have 

U = 1 2 P~t + U , (3.7) 

which is the light-cone gauge. In that case, the constraint Too = gives, upon substitution of the 
momenta, 

- 2P + P_ = PlAijXtXi + + d ° Xl d° X \ (3. 8 ) 

V is given implicitly by the constraint, leaving the transverse coordinates as the physical degrees 
of freedom. Taking the metric function Aij = XujSij, and defining the dimensionless mass Mu\ = 
H(i}l 2 P-, we have 



□l i T M|r = {\ {i) =±fa). (3.9) 



Consider, first, the zero modes: 



a' 



Xlir) = x l cos {M (i) r) + — - sin (M {i) r) p j (A w = fa), (3.10) 

W 

X*(t) = x l cosh (M (l) r) + -?L sinh (M (i) r) p* (A (i) = -^)), (3.11) 

(») 

which of course match the geodesies (2.14) and ( 2.15j ) above. The centre-of-mass of the string 
will thus oscillate about the origin of those directions associated with positive eigenvalues, and will 
accelerate forever in the directions associated with negative eigenvalues. 
The non-zero modes are solved for by making the Ansatz 

X\r,a) = Y^Cy {uJnT+na \ (3.12) 
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giving 

c4 = n 2 ±Mg } (A (i) =± M J)), (3.13) 

For negative eigenvalues, the set of modes with |n| < M(j) have lo 2 < 0. As in the more general case 
considered in |28|| , some of the uj n become imaginary, and would seem to correspond to unstable 
modes. However, this is an instability with respect to mode creation on the string worldsheet, in 
precisely the sense of ^ |||, and not in the sense that we are concerned with here. We do not 
expect to see truly stringy instabilities in our perturbative analysis. Again as in j^], these stringy 
instabilities exist only for some range of n. In particular, we have unstable modes only if 

fi(i)P- > T, (3.14) 

T = 1/lg being the string tension. Unstable modes only appear if either the light -cone momentum 
of the string, or the curvature of the background, is large with respect to the string scale (the same 
is true in the more general case considered in p8[ ). 

In the directions associated with positive eigenvalues, the oscillator modes are given by [10] 

Xi = E — « ema + <e- irur ) e~^, (3.15) 



where, for n > 0, u> n = x /n 2 + M?~ and, to ensure reality, 



|t e -W\ e n n r 



cj_„ = -u n , (a l n ) = a l _ n , (a % n ) = a l _ n . (3.16) 
In the remaining directions, associated with negative eigenvalues, we have 

xi = n/r E J- « eina + ^~ ina ) ^ nT + M E fr + 

\n\>M {i) n \n\<M (l) " 

(3-17) 

where, for n > 0, io n = ^Jn 2 — M 2 ^ , and the same conditions as in ( |3.16|) apply to ensure reality. 
The (3 l n operators are associated with the exponentially growing/decaying modes and we have set 
uj n = i£l n for | tt. | < M(j) where, for n > 0, il n = M 2 ^ — n 2 . Then the expansion for the unstable 
modes is manifestly real. To include all the relevant modes in this case, we also set f2_ n = — Q n . 
The non-vanishing Poisson brackets are then 

[x\pi] pB = 5 i i, 

[ a lni a n]PB = [a % m ,a? n ]pB = -WmW.O^, 
\? m ,Pk]PB = [(/4) f , (/%) + ]PB = "mWn.oP. (3-18) 

Upon quantisation, the modes with real co n will give rise to the usual harmonic oscillators, but 
those with imaginary iv n have an extra factor of i. In other words, the fi l n do not correspond to 
harmonic oscillators, rather they are associated with the unstable worldsheet modes. 
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Substituting for the mode expansions of the worldsheet scalars and the respective momenta into 
|) gives 

-2(2vr) 2 P + P_ = (2^) 2 P 2 A (j) xV + pV + A E + 



a' 



+ ^ E Ot^ + M)' (£•»)*)> ( 3 - 19 ) 



a 

|n|<M {i) 

where the first oscillator sum runs over all directions associated with positive eigenvalues, and 
over those modes for which |n| > Mu\ in the directions associated with negative eigenvalues. The 
second sum applies to the latter directions only. The P 2 contribution comes from the zero modes, 
so would also be present in the case of a particle. In either case, it is not at all obvious as to how 
this dispersion relation gives us information on what our flat -space intuition would call a "mass" . 
In flat space, we have 

- 2P+ = — (pV + oscillators) , (3.20) 

which gives well-defined behaviour; in particular, P + — > as P_ — > oo, so that we never have 
faster-than-light propagation. In the case at hand, however, we have 

- 2P+ = P_A (i) zV + — (pV + oscillators) , (3.21) 

the non-standard dependence on P_ giving rise to "unstable" behaviour: depending on the sign 
of A(j), we have P+ — ► ±oo as either P_ — ► oo or fj,^ — > oo. In either case, we know that in this 
limit, only the zero modes survive, since all the stringy modes with n > are unstable; the string 
is literally ripped apart. Note that modes with small P_ behave as if they have a positive mass, so 
in this infra-red regime, there are no tachyonic excitations. 

At any rate, based on flat-space intuition it would appear that we either have slower-than- 
light or faster-than-light propagation depending on this sign. Negative eigenvalues seem to give 
tachyonic behaviour, but this is presumably to be expected since particles (and strings) accelerate 
forever. Even so, we will not find any unstable modes in the analysis below. As we will see, the 
problem in this simple case reduces to the quantum mechanics of the inverted oscillator in which 
similar behaviour is observed: quantum particles accelerate forever, but the system is perturbatively 
stable. Various subtleties concerning this dispersion relation are addressed in , where it is argued 
that imposing boundary conditions at some finite x t x t = I? effectively removes the problems 
associated with large P_. 



To make contact with (1.1), we rescale 

< = v^4 aLn = V^ai, (3 l n = 0L n = (3.22) 

and likewise for the a and j3\ and combine the zero modes as 

° = 7S^ ( " y - iM <->^ ^ = ^kp w + ,M{i,x ' y (3 ' 23) 
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Upon quantisation, the new operators obey 

[4,4) = ^, [ai n ,ai) = 8 mn S l \ [b i m ,V n )=i6 mn 5 i i, (3.24) 
and likewise for the a and f3' , and we have 

-2P + P- = ±J2[ M W N o ) +Y,^ )N % )+ E ^W], (3-25) 

a i \ «>0 \n\<M (i) J 



where M {i) = ±/i (i )/ 2 P-, < = ^n 2 + M 2 } , AT® = <4 + atf n and = b\b\ + (b^t^Jt. 

4 Linearized field equations 

We will consider perturbations of the type IIB version of the vacuum plane wave discussed in section 
2. We have some results regarding the non-vacuum solution, but will not discuss them here since 
this case is much harder to analyse in any detail. Moreover, the vacuum solution is the cleanest in 
which possible negative (mass) 2 terms are generated on the string worldsheet. 

The field content of the bosonic sector of type IIB supergravity is: a complex scalar, $> = 
X + ie~^, combining the axion and dilaton, % and (/>; a complex 2-form potential with field strength 

G 3 = dA 2 = H 3 + iF 3 , (4.1) 

combining the NS-NS and R-R 3-form field strengths, H 3 and F 3 ; and an R-R 4-form potential 
with self-dual field strength 

F 5 = *F 5 = dC A - 1q (A 2 A G* 3 ) , (4.2) 

o 

where * denotes complex conjugation and * denotes ten-dimensional Hodge duality. We take the 
field equations from [61], perturbations of which, in a vacuum background, simplify dramatically. 



4.1 Scalar perturbations 

Throughout the analysis we will find ourselves solving the field equations similar to those for a 
massless scalar, so consider the complex scalar, which satisfies 

□$ = 0. (4.3) 

We want to look for instabilities, i.e. modes which grow exponentially in time, where "time" is u. 
We will also demand that the modes are normalizable — either L 2 normalizable or delta-function 
normalizable — in v and all transverse directions. We have already discussed above that we can 
always decompose the modes with respect to the Killing vector d/dv, so we make the Ansatz 

$ = e i( - p - v - p + u ^(x), (4.4) 
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where J-{x) = 4>i(x) + i(f>2(x) depends on all transverse coordinates. Then 

□* = e i{ - p - v - p+u) [(2p_p+ + H(u, x)pl)F(x) + d k d k T{x)} . (4.5) 

We assume that p- is real, else the modes will not be normalizable in the v direction, and unstable 
modes will have imaginary p + . Our approach will be to consider those modes with real p+ and 
show that in fact this includes all modes with real momentum p_ , so 

[2 P -p+ + H(u, x)p 2 _)<t>{x) + d k d k ^{x) = 0, (4.6) 

for both <f> = <pi and <fi = fa. We will see that precisely the same equation controls the perturbations 
of all supergravity fields, and will analyse the solutions in the following section. 

4.2 Metric perturbations 

As usual in perturbation theory, we need to choose a gauge to restrict ourselves to physical degrees of 
freedom. Although the standard transverse, trace-free gauge considerably simplifies the linearized 
field equations, this choice does not entirely fix the gauge. We will rather make use of the light-cone 
gauge, which is well-suited to our background spacetime; but also has the advantage of fixing the 
gauge freedom entirely. The field equation^] is 5R a i, = 0. 
The light-cone gauge is defined by 

h va = 0, V a, (4.7) 

and substituting for this gives 

- 25R UU = Uh uu - 2d u (d v h uu + dih iu ) - didjHhij - 2d i Hdjh ij , (4.8) 

-25R UV = -d v (d v h uu + dihi u ), (4.9) 

-25R ui = Dhui + djHdvhij -di(d v h uu + djh ju ), (4.10) 

-25R VV = 8$ha, (4.11) 

-25R vi = -d v (d v hi u + djhij), (4.12) 

-25Rij = Dhij. (4.13) 



For the maximally supersymmetric plane wave, these results were effectively derived in [10, 33 1. As 
in that case, the equations 5R va = appear as constraints, which we solve by taking 

ha = 0, d v hiu + djhij = 0, d v h uu + dihiu = 0, (4.14) 



7 One often thinks of the linearized vacuum field equation to be the vanishing of the Lichnerowicz operator [K| 
acting on the metric perturbations. This is true only in the transverse, trace-free gauge however. That is \7hh a b = 
— 2SR a b only if h — g ab h a b = ^bh h a = 0. Since we will not be working in this gauge, our vacuum field equation is 
SRab = 0, and not = 0. 
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which leaves the equations 

Dh tj = 0, (4.15) 

Dh iu + djHd v hij = 0, (4.16) 

Uh uu - didjHhij - 2d i Hd j h ij = 0. (4.17) 

We thus need only solve 

Uhij = 0, (4.18) 
all other components of h a b following from via the constraints ( [4.14 ). Working in the light-cone 



gauge reduces the degrees of freedom to the physical ones only, the transverse components hij, 
minus the trace. 

We make the Ansatz 

hij = & {tij e^- v -P^F(x) } , (4.19) 

where, as in the previous subsection, T(x) = 4>\{x) + ifaix). £y = ^nj\ denotes a constant 
polarisation tensor which satisfies £a = 0. Again we consider modes normalizable in v (real p_) 
with real p+. Then 

Uh^ = K {^e^-^ [(2p^ P+ + H{u, x)p 2 _)F{x) + d k d k H*)] } = 0. (4.20) 
Writing this in terms of real and imaginary parts, we have 

[2p_ P+ + H(u, x)pl]<P(x) + d k d k <p(x) = 0, (4.21) 



for both cj) = 4>\ and <p = 4>2, just as in ( |4.6| ) above. 
The constraints ( 4.14| ) give 



h iu = ^l^Cij e^- v -^djT{x)^ , h uu = e^-v-^didjH*)} , (4-22) 

though we still need to check for consistency with the field equations ( |4.16| ) and ( 4.17| ). Substituting 
for the above results, we find 

(4.16) 5ft e i( P-v-P+u)dj [(2p_p+ + H(u,x)p 2 _ )T(x) + d k d k F(x)] 1=0, (4.23) 

(4.17) R j^-fcj jto-v-P+^didj [(2p_ P+ + H(u, x)p 2 _)F(x) + d k d k T{x)] 1=0, (4.24) 



so consistency is guaranteed and we need only take 4>\ and 4>2 to satisfy (4J3). Solutions of this 
equation thus fix all metric perturbations. 
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4.3 4— form perturbations 

In addition to ( |4.7| ), we now also impose 

SC vabc = 0, V a, b, c, (4.25) 

so that 

8F vai .. M4 = d v SC ai ... a4 . (4.26) 
Perturbing the self-duality condition F§ = -kF§ gives 

8Fvi 1 ...i 4 = — Tj £ ii...i4ji— i4^uii— j4) (4.27) 

&F uv ij k — ~^Sijkii...i^^Fi 1 (4.28) 
8F u i x j 4 = ^y£i 1 ...i 4 j 1 ...j 4 5i ; i t jf 1 ...j4 , (4.29) 



where e^...^ = ..j g - The 4-form potential is anti-self-dual in eight dimensions, since ( 4.27 ) 
gives 

^Cil...l4 = ~-^^ii...i4ji...j4^Cji...j4- (4.30) 



Furthermore, ( [4.28 ) is just the constraint 

d v SC uij k + didCujk = 0, (4.31) 

which can be solved for SC u ij k . 

The 4-form equation of motion is 

V e 5F eabcd = 0, (4.32) 
the {u, v,i, j} component of which can be solved by taking 

d k SC kuij = 0, (4.33) 



and the {v,i,j, k} component gives the constraint (4.31). The remaining components are 

D5C ijkl = 0, (4.34) 

U5C uijk + mdvdCujk = 0, (4.35) 

and, again, we will see below that the solution for 5C u ij k found from the constraint ( |4.31| ) is 
consistent with this field equation. 



Since the 4-form is real, we make a similar Ansatz to ( |4.19 ) above: 



SC ijkl = R {eyw e i{ P- v ~ p + u ^{x)] , (4.36) 
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where 

1 

eii...j 4 — — Tj^ii—HJl— H e 3l— 34,1 

is a constant polarisation tensor, anti-self-dual in eight dimensions, and F{x) = 4>i(x) + i(j>2{x). 
Then 

D8C ijkl = ft {e^ e *(P-«-*M-«) [(2p_p+ + x)p 2 _)^(x) + d m d m T(x)} } = 0. (4.38) 
The constraint Q4.3l|) gives 

SC uijkl = &{j-e lijk e^-^+^i^)} , (4.39) 



so that ( 4.33; ) is an identity, and substituting this into the field equation (4.35) gives 



^\p~ eiijk j {p - v - p+u) di [(2P~P+ + H(u, x)p 2 _)T(x) + d m d m T(x)] |=0. (4.40) 
Again, both <f> = (pi and (j) = 4>2 must satisfy ([4 . 6|) . 



4.4 2-form perturbations 

The 2-form perturbations go through in the same way as the for the 4-form above, the only 
difference being that they are complex. Taking 

5A va = 0, V a (4.41) 

gives 

5G vab = 8 v 5A ab . (4.42) 

The equation of motion is 

V c 5G cab = 0, (4.43) 
the {u, v} and {v, i} components of which are constraints which can be solved by taking 

diSA iu = 0, d v 5A ui + djSAji = 0. (4.44) 

This leaves 

DSAij = 0, (4.45) 
DSAui + djHdJAji = 0, (4.46) 
the structure of which is the same as for the 4-form perturbations. Writing 

6 Ay = Cij e^- v ~P+ u ^(x), (4.47) 

for some constant complex polarisation tensor Qj = C[ijh gives 

DSAij = ( tj e i (P—P+«) [(2p^ P+ + H(u, x)pt)F(x) + d k d k F(x)] = 0. (4.48) 
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The constraint (4.44) gives 



SAui = —Cji Jto-v-P+^djFix), (4.49) 



and the field equation ( 4.46 ) is 



e i<P-«-P + ») dj [( 2 p_ P+ + H(u, x)p 2 _)F(x) + d k d k F(x)] = 0, (4.50) 
so we see precisely the same structure as above. 

5 Classical stability 

In the above section, we showed that the linearized field equations effectively reduce to the problem 
of solving (|4.6| ): 

[2p_p+ + H(u, x)p 2 _]<t>(x) + d k d k <j>{x) = 0. (5.1) 

Solutions of this equation control the perturbations of all supergravity fields. So far we have 
considered those modes with real momentum p_ that also have real p + . Here we will see that in 
fact all modes of real p_ have real p + . This is done by finding the general solution to ( |4.6| ) which, 
using a separation of variables, comes down to a second order o.d.e. with two known, linearly 
independent solutions. Note that we have ignored an arbitrary function of u which could be added 
to the Ansatze above. That is, we could have taken 

$ = e i( P- v - p+u ^(x) + /(u), (5.2) 

and likewise for the other supergravity fields. This is a manifestation of the fact that u = constant 
is not a Cauchy surface, but a null surface. However, treating it as a substitute Cauchy surface only 
misses those modes which travel along constant u |5^| , parallel to the wave front. The presence 
of such a function f(u) is familiar in light-front studies |63| and is determined by the boundary 
condition that fields should vanish asymptotically, requiring f(u) = 0. 



Let us consider the simplest case, taking Aij = \(i)dij with the eigenvalues = /if-j, 



1, . . . , n and A(j) = -/iLi = n + 1, . . . , D — 2. Separating variables as 

Fix) = F 1 (x 1 ) . . . F^Giiy 1 ) . . . G D ^ n (y D - 2 - n ), (5.3) 
where we denote the directions associated with negative eigenvalues by {y l , . . . ,y D ~ 2 ~ n }, gives 

F(i) + " V 2 {l )P-(x 1 ) 2 ) F {t) = 0, G'fa + (E {i) + 4 /_(yT) % } = 0, (5.4) 



where 

D-2 



^2 E {i) =P+P~- (5-5) 
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The equation for the Fm is just the Schrodinger equation for a simple harmonic oscillator, whereas 
the equation for the Gm has an inverted oscillator potential. Of course, this is the origin of the 
oscillatory and exponential behaviour in the geodesies and zero modes of the string as discussed 
above. 

Asymptotically, as \x\, \y\ — > oo, we have 

F{{){x i) ^ e ±^ )P -(^)\ Q^yi) ^ _L cos 1 gin (l^p^f^ , 

(5.6) 

for each F^(x % ) and G^{y % ). In the x 1 directions (with positive eigenvalues), the modes are 
either exponentially damped, or are exponentially growing — these are the normalizable and non- 
normalizable modes of |35[ |. In the y l directions (with negative eigenvalues), however, the solutions 
are oscillatory, the two signs giving left- and right-moving plane waves. The former are easy to 
deal with: taking 

F ( , ) (x i ) = e ± ^P-^ 2 f (i) (x i ), (5.7) 

gives 

/g } ± 2^ (i) p_x i /( i) + (E (i) ± = 0, (5.8) 

which is of course solved in terms of Hermite polynomials, up to an overall constant, giving a 
normalizable solution [f3^, 33] 



F^x*) = c n e-^v^) 2 H n {^jr^Zx l ), n=\ ( - l) . (5.9) 

z \H i (i)P— / 

For real n, i.e., for real Eu\, these solutions furnish a complete basis of normalizable modes solving 
the equation of motion in the x % directions. 

We now turn to consider the directions associated with negative eigenvalues. Upon defining 
z % = y i /a^, where a%, = 1/(2/^) \p-\), the equation for the G^ in ( |5.4| ) becomes 

1 



+ ( 4^) 1^(0=0, (5.10) 

where = —E^c&y This can be solved in terms of parabolic cylinder functions giving, for each 
z l , two independent solutions |34|] 

Mz)=W(a,z), Xa (z) = W(a,-z). (5.11) 

The first solution has the limiting behaviour 

l~2k 

ip a (z -> oo) -> y — cos(z 2 /4 - alogz + 7r/4 + v?/2), (5.12) 



^ a (z^-oo) -> y^|-sin(z 2 /4-alog|z| +7T/4 + V3/2), (5.13) 
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where ip = argr(^ + ia) and k = yl + exp(27ra) — exp(7ra), and the second behaves as 

T 




-oo 



sin(z /4 - a log z + 7r/4 + <p/2), 



/2A; 

— sin(z 2 /4 - a log 1 2 1 + tt/4 + ip/2), 



To show that these functions are orthogonal, we write ( |5.10 ) as 

G'a + ~^z 2 G a = aG a , 



(5.14) 
(5.15) 

(5.16) 



where G a denotes either of the two functions in ( 5.11 ), multiply both sides by Gb and integrate by 
parts, giving 



G a (z)G b (z) dz = 
l {a - b) 



\GbG' a — GaGfj 



b\-L 



Taking the large L limit and using the limiting behaviour discussed above, we find 

/oo 
Xa {z)Xb{z) dz = 2irVl + e 2 ™ 5 (a - b). 
-oo 



(5.17) 



(5.18) 



Substituting instead for G a = ip a and Gb = Xb gives 



1 



ipa(z)xb(z) dz = 
L (o - b) 



cos [(a — b) log L]. 



For a = b + e, with e small, 



which shows that 



k a = h + ire 



„2tt6 



Vl + e 2 ^ 6 

^a(^)Xfc(^) 



7re 



?ab, 



Vl + e 27ra 

where <5 a fc is the Kronecker delta. We can then write an arbitrary function as 



G(z) = J da[a{a)il) a (z) + (3(a)x a {z)} 



where a and (3 are found using the orthogonality relations 
1 



a(a) 



2tx^J\ + exp(27ra) 



dzVa(^)G'(z) 



/3(a) 



(5.19) 

(5.20) 
(5.21) 

(5.22) 

dz Xa (z)G(z). (5.23) 



27ryl + exp(27ra) 

To recap, in the directions associated with negative eigenvalues we found solutions in terms of 
cylinder functions and in the directions associated with positive eigenvalues we found solutions in 
terms of Hermite polynomials. If the parameters a and n are real then these form a complete basis 
and we may use the orthogonality of these functions to decompose any initial data set in terms of 
modes with real a, n. However, real a, n implies real Eu\ and with it, from Q5.5| ), real p+. That 
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is, an arbitrary initial perturbation profile is composed of modes with real p+, and so the system 
would seem to be stable. 

To clarify the above, let us consider what happens in the more familiar case of a free scalar field 
with negative (mass) 2 : 

- 8$<p + d 2 x <p + d 2 <p + raV = 0. (5.24) 
To follow an analogue of the above argument we consider a mode of momentum p x , 



cp(t, x, y) = A(y) cos(ujt - p x x) + B(y) cos(wt - p x x), (5.25) 

giving 

d 2 A(y) = -(oj 2 +m 2 - p 2 x )A(y) = -p 2 y A(y), (5.26) 

with a similar equation for B(y). The solutions to this are either exponential or trigonometric, 
depending on the reality of p y - for normalizable solutions we need to take real momentum, p y , 
giving us trigonometric solutions. We are now able, using standard Fourier analysis, to decompose 
an arbitrary initial perturbation in terms of this basis of real p x , p y . This is the analogue of being 
able to write solutions to ( |4.6[ ) in terms of cylinder functions with real o and Hermite polynomials 
with real n. The difference now is that uj 2 = p 2 — m 2 , so small momentum (long wavelength) states 
have imaginary uj whereas for the plane-wave we found that p + was always real. 

One should be careful, however, since there are subtleties associated with the lack of a true 
Cauchy surface (some of which are discussed in although the same conclusions are reached 

in that work. As we have mentioned above, it would seem that evolving initial data from a 
u = constant surface only misses modes which are unlikely to exhibit unstable behaviour, since 
they are unaffected by the wave. 

6 Conclusions 

The propagation of a closed string in the light-cone gauge on the vacuum plane-wave gives rise 
to world-sheet scalars with a tachyonic mass. In particular, the zero mode of the string, which 
determines the supergravity modes, suffers from such a mass term; this motivated the study of 
the classical stability of the supergravity solution. Although such spacetimes are super symmetric, 
and supersymmetry is usually understood to imply stability, the standard arguments do not apply 
as the metrics are not asymptotically flat (they are Lorentzian symmetric spaces). It may turn 
out however that such methods could be modified to cover the plane-wave, and this would be an 
interesting avenue of research. 

Taking the IIB equations of motion, perturbed around the vacuum plane wave, we have shown 
that in a gauge suited to the spacetime (light-cone gauge) the system boils down to the study of 
a single equation, (4.6). The simplicity of this equation allowed us to find its general solution in 
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terms of known functions. Moreover, an arbitrary initial perturbation was seen to be decomposable 
in terms of modes with real p+, corresponding to oscillatory motion, and thereby providing strong 
evidence for the stability of the spacetime. 

Of course, one would like to be able to prove that instabilities do not exist, by showing that 
there are no normalizable solutions of ( |4.6| ) with imaginary p + . Unfortunately, however, we have 
not been able to do this. 

On general grounds, one might suspect that the presence of tachyonic worldsheet scalars would 
not give rise to any instabilities of the background, for the simple reason that the mass of these 



scalars is purely an artifact of choosing the light-cone gauge on the worldsheet [25, 39 1. Choosing 
instead, for example, the static gauge, gives rise to a worldsheet theory of massless but interacting 
scalars [^9[. Indeed, simply changing from Brinkman to Rosen coordinates in the background 
spacetime has a similar effect on the worldsheet. It would be interesting to investigate how the 
interactions in different coordinates or gauges mimic the effect of a mass, tachyonic or otherwise. 
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